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Abstract. Many cosmological observations call for the existence of dark matter. The most 
direct evidence for dark matter is inferred from the measured flatness of galactic rotation 
curves. The latter is based on Newtonian gravity. Alternative approaches to the rotation 
curve problem by means of general relativity have recently been put forward. The class of 
models of interest is a subset of the axially symmetric and stationary solutions of Einstein's 
equations with rotating dust. As a step toward the understanding of general relativistic 
galaxy models, we analyse rigidly as well as non-rigidly rotating (Post-)Newtonian space- 
times. We find that the Newtonian limit of the considered general relativistic galaxy model 
leads to Post-Newtonian terms in the metric. 



1. INTRODUCTION AND MOTIVATION 

Miscellaneous physical observations support the postulation of a new component of 
matter, non-luminous and only interacting through gravity. The phenomena accord- 
ing to which dark matter is postulated manifest themselves on multiple cosmological 
scales, ranging from galaxies up to the primordial Universe. The concordance model 
that enfolds observations of the CMB as well as of the large scale structure currently 
charges a value of Sldm — 0.23 (WMAP 2008). Herculean efforts are being made in 
order to find an explanation for the dark matter findings. The attempts range from 
modifications of the Einstcinian, and therewith of Newtonian, gravity to extensions 
of the standard model of particle physics that involve new, yet undetected particles 
that could do the job. But it is important to stress that there is no a priori reason 
to believe that all of the missing matter problems on many different physical scales 
have a common explanation. 

Here we will be concerned with the dark matter problem on galactic scales. The 
missing of a Keplerian velocity fall-off (w(r) oc l/v^) in most of the observed rota- 
tion curves represents the main classical evidence for the existence of dark matter in 
galaxies (Rubin et al. 1980, Binney & Tremaine 1994). There appears to be far more 
matter in a galaxy than the luminous part suggests. From the flatness of the rotation 
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curves up to maximally measurable radii, the existence of vast halos of dark matter 
is concluded. This leads to a cosmological lower bound of fldm ^ O.l (Amsler et al. 
2008). The famous missmatch occurs when comparing the observed rotation curves 
to Newtonian gravity. But the Newtonian limit appears well justified for a galactic 
system since the gravitational field strength and the occurring velocities are small. 

However, in a recent (Cooperstock & Tieu 2005-2007) series of works Cooperstock 
and Tieu (CT) suggest a new approach to the galactic dark matter problem, namely 
through general relativity (GR). According to these authors the application of full GR 
to galaxies introduces new effects, cumulating in the bold statement that dark matter 
can be completely obviated. What may be the motivation for such a try if velocities 
and fields are small? According to CT, already in a stationary system that is purely 
gravitationally bound, inherent non-linear terms occur in GR whose magnitudes may 
not be negligible. Unlike in the Solar System where GR effects are important but 
minute, CT argue that the constituents of a galaxy, the stars, are not test particles in 
the field of a central mass but both constitute and follow their mutual common field. 
This relation does of course not increase the field strength but it shall influence the 
dynamics of the system in the proper theory. Tangibly, CT fitted flat rotation curves 
in the GR model and are able to return realistic density profiles using luminous matter 
only. The resulting integrated masses of CT galaxies lie below the Newtonian values 
but are larger than the total masses derived from Modified Newtonian Gravit}0. 

Numerous and severe criticism has been published against the concrete CT model 
in the meantime. In several works it has been pointed out that there exist patholo- 
gies within the energy momentum tensor or various other imponderabilities with the 
model. For instance Korzyhski argues that (A) a proper fiat asymptotical limit is not 
contained in the CT model because of the considered gauge; and that (B) the CT 
model must be unrealistic because unexpected additional matter sources at z = can 
be found (Korzyhski 2005): a singular behaviour of the energy-momentum tensor in 
the plane is revealed because of the existence of a residual non- vanishing Komar mass 
within a volume that has been shrunk to zero. 

Even though the concrete model of CT is problematic the main idea points our 
attention to GR modelling of galaxies which is a valuable point. This is an important, 
yet difficult field where less progress has been made than one would expect. There 
exist only a handful of known axially symmetric solutions in GR and all of them are 
not viable as a galaxy model, see e.g. (Islam 1985). The only existing global solution 
that comes near to this is the one by Neugebauer and Meinel for a disc of rigidly 
rotating dust (Neugebauer & Meinel 1994, 1995). This solution was found by the 
inverse scattering method and its formulation is exceedingly complicated and mostly 
relevant for numerical modelling. Instead of searching for an exact solution resembling 
a galaxy, we study the Newtonian limit and the Post-Newtonian contributions to 
axially symmetric and stationary spacetimes. 

2. AXIALLY SYMMETRIC AND STATIONARY SPACETIMES 

The most general axisymmetric {dipg^^i, — 0) and stationary {dtg^^v — 0) spacetime in 
GR takes the form (Stephani et al. 2003) 

ds^ = [gMNdx^'Ax'' -f W^dip^] - e^"" (dt + Adipf , (1) 

^For original refs. see (Milgrom 1983) and for the relativistic extension (Bekenstein 2004). 
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with notation {x'^ , , x'^ , x^) = {t,r,(p,z) and M,N take values 1 or 3 . It is always 
possible to go to isotropic coordinates gMN = c^^^mn without touching the generality 
of ([!]), c.f. (Stephani et al. 2003). So we get the final form of the most general 
axisymmetric and stationary spacetime named after Lewis and Papapetrou (LP) 



ds2 = e-2C^ [e2fc(clr2 + dz^) + W^d^""] - e''' {dt + Ad^f 



(2) 



with the free metric functions J7, fc, W, A being solely functions of r and z . 

It is possible to simplify ^ a bit more, but only under certain assumptions. We 
will demonstrate that if and only if the metric function W is harmonic it can be 
transformed tcH W — r . Consider a complex coordinate transformation f{r + iz) = 
W + iV introducing an additional potential V . Then we have from p = W{r, z) and 
h = V{r, z) the differentials dp — 



or 



dW 
dz 



dz and dh = ^dr 



az 



With the 

coordinates p, h being only dummies, introduced for bookkeeping reasons. We insert 
the transformation into , written in terms of p, h , and have 



ds 



e'^idp' + dh^) + p'dip^ - e^^ {dt + Adip) 



dVdV , , 
2— — drdz 
or az 



\ dr J \ dz 



dz^ + 2— — drd2 

or az 

+ W^dip^ \ - e^^{dt + Adipf . 



dV 
dr 



ds' 

2 



dr^ 




(3) 



Requiring formal invariance as compared to the original metric, we see that the mixing 
terms should vanish. That is exactly provided by the Cauchy-Riemann equations for 
W and V 

dW dV , dW dV 

-7r = ^ ^"^^ -ir = -IT ■ 

or oz az or 

Moreover, with the help of the Cauchy-Riemann equations, we see that the coefficients 
of dr^ and dz^ can be combined to a positive definite quantity 



dr ) 



V dr , 

=2fc i 



' - (dWV 



. dz ) 



= K>0, 



(5) 



such that we can combine e^'^K = e'^ and so obtain ^ via Q. Thus we have 
shown that it is possible to simplify the general LP form ^ by allowing W = r , 
which is only possible if the transformation / is analytic, that is W (and also V) 
must be a harmonic function, A'^'^W = . This condition for W holds for exterior 
(vacuum) solutions that are stationary and axisymmetric (Islam 1985, Stephani et 
al. 2003), thus any GR model of a galaxy must approach this class of solutions where 
the matter density drops to zero. Then, we can write down the LP metric in isotropic 
coordinates (or Weyl gauge) 



-2(7 l„2k 



{dr' + dz') + r'dip' 



o2C/ 



(dt + Adipf 



(6) 



Also W = 1 is possible then, but this case is of no interest in our case. 
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The spacetime applied in the CT model is the following 

ds2 = e''-™ (dr^ + dz^) + r^e-'^'dif'^ - {dt + Nd^f , (7) 

together with corotating dust as the matter. Obviously, the CT metric does not 
belong to the class of the most general stationary and axisymmetric spacetimes; it 
belongs to the subclass of LP solutions in the Weyl gauge, and it is not obvious that 
this is justified in the presence of dust. 

Now we ask what solution could potentially be a Newtonian counterpart to the 
CT model ([7]). The 'Newtonian approximation', that is the metric that reproduces 
Newtonian physics, is given by, see e.g. (Misner, Thorne & Wheeler 1973), 

ds^ = -(1 + 2(t>)dt^ + dr^ + r^dif^ + dz^ , (8) 

where <j){r, z) is the Newtonian gravitational potential. CT work in a comoving frame, 
thus we consider the rotation of the Newtonian metric. For simplicity, we start with 
rigid rotation in ([5]) via ip — Lp' — Lot. The exact result can be brought to the form 



(1 + 20-^2^.2) ^ 



ds^ = (dr^+dz^)+ ^l'^'^'^, ,, rW-(l+20-c^V) 

(9) 

This is the rigidly rotated Newton metric. In this form we can directly compare it 
with the LP metric in Weyl gauge ®, and we notice a discrepancy at linear order 
in (f> , looking at the dip^ term. Interestingly, the rigidly rotated Newton metric ^ 
is not consistent with the Weyl subclass of the LP solution ©. Note that ^ is in 
perfect accordance with the general form of the latter 

One could now ask whether the situation might be easily cured with the help 
of a coordinate transformation. We show that this is not possible. Above we have 
derived the exact conditions under which the general and the isotropic LP metric can 
be transformed into each other: the function W must be a harmonic function with 
respect to the two-dimensional Laplacian A^^)]^^ = q. In the present case, of the 
rigidly rotating Newton metric, W is given by 



W ^r^/l + 2(j). (10) 
Expanding to linear order and applying the Laplacian yields 

A(2)M^ = rA(3)0 + </>,^ =47rGpr + 0,r- (H) 

Note that we can use the Poisson equation because the potential is Newtonian. After 
repeating some facts from potential theory we will show that A^^)^/ in fact does not 
vanish in general. 

Given the general problem of solving the Laplace equation with the appropriate 
boundary conditions for a disc-like distribution of matter, the solution for the potential 
can be obtained via separation of variables, c.f. (Binney & Tremaine 1994), 

roc 

(j){r,z)^ S{k)Jo{kr)e-''^'^dk. (12) 
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A given surface mass density S(r) is then characterised by the according Hankel 
transform 

/•oo 

S'(fc) = -27rG/ Jo(fcr)E(r)rdr. (13) 

Now we can use these expressions for the evaluation of (fTT|) . 

Case (A) z — Outside the disc the Newtonian potential fulfils the Laplace 
equation, such that the expression (jlip takes the form 

/•oo 

A(2)w^ = -/ 5(fc)Ji(fcr)fce-'=l^ldfc at z^Q, (14) 
Jo 

which will not vanish in general. As a simple example we consider the Mestel disc 
model of a galaxy (Mestel 1963). In the Newtonian Mestel model a flat rotation curve 
can be reachecH, although for the price of an infinite total mass. The Mestel model 
is characterised by a surface mass density that falls off inversely with the distance 

I](r) = ^. (15) 

In a Mestel galaxy the surface density Hankel-transforms as S{k) = —2TTGY,oro/k. 
Using this we can integrate directly and obtain 

A^^^W = 2TTGJ:oro(-~^===] at z^O. (16) 

Case (B) z ~Q — We want to show that (fTT|) is non-zero also here. Let us assume 
the contrary and see what happens. If we assume that H^S^^W = was true then 
equation (|lip gives an identity. This we integrate over z for some e > and then 
revoke the operation by performing the appropriate limit 

/£ PS POO 

(5(z)S(r)dz = lim / / S'(fc)Ji(A:r)fce-'=l^ldfcdz. (17) 

Because the exponential term on the right hand side acts like a damping factor, the 
modulus of the integrand will reach its maximum at z = 0. Thus, as an upper 
estimate, we can set the integrand of the right hand side to be constant in z and 
therefore the integration and limit procedure give zero. Then, for all other z the 
expression will be zero more than ever and we obtain 47rGS(r)r ~ 0. This will not 
hold generally for any realistic model, hence producing a contradiction, and therefore 
A'^^^W{r, z) = is not true at the surface z = either. 

Interestingly, the pure Newton metric ([5]) cannot be made compatible with the 
LP metric in Weyl form and thus is also not compatible with the CT model. As all 
exact vacuum solutions fall into the Weyl class as well, it seems that the Newtonian 
metric is not appropriate to describe the physics of rotating systems at leading order. 
We should go one step further and consider the Post-Newtonian (PN) metric 

ds^ = -(1 + 20)di2 + (1 - 2?A)(dr2 + r^d^^ ^ ^^^2^ ^ ^^g^ 



■^The flat rotation curve in the Mestel model is obtained from the Hankcl transform of II15I I. 
inserted into the formula for the rotation curve: v'^{r)MeB = r{d4>/dr)z^Q = 27rGSoro . 
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with some additional PN potential ^ . Sometimes, this metric is referred to as the 
'longitudinal Newton' metric. The reason might be that the order of magnitude of the 
coefficient of the spatial part and the order of the Newtonian correction are the same. 
Nevertheless, conceptually this makes an enormous difference. In classical Newton 
Gravity there exists no curvature of space, the three-space is always euclidian. This 
is exactly reflected in the Newton metric ^ and therefore we refer to (fTS]) as the PN 
approximation; for an extensive discussion see e.g. (Misncr, Thorne & Wheeler 1973). 
Now, applying stiff rotation to the PN metric results in (including all orders) 

(1 + 20)(1-2V;) 



ds^ = (1 -2V')(dr^ +dz-') + 



[1 + 20 - (1 - 2V')wV2 



[1 + 20 - 



dt 



(1 - 2V')w2r2 

^_(l^-2V0wr^__ 

[1 + 20 - (1 - 27/>)w2^2 



■ dip 



(19) 



For the simplest PN case, i.e. = ■0 , this is consistent - neglecting terms 0(0^) 
with the Weyl gauge © of the LP class via 



e^" = (1-20)[1 



20- (1 
(1- 



- 20)^2, 
20)wr2 



e^^ = [1 + 20 - (1 - 20)wV^] 



[1 + 20 - (1 - 20)w2r2] 



(20) 



Thus the PN metric after rigid rotation belongs to the class of isotropic (Weyl) ax- 
isymmetric and stationary solutions, if = -0, whereas the rotated Newton metric 
^ does not allow for that simplification. 

Differential rotation is more realistic. So next we relax the condition of rigid 
rotation and consider transformations — if' — uj{r)t with the the PN metric, giving 



ds^ 



(1 - 20>)dz2 + (1 - 2i!)(l + r^ijlf)dr^ + 



(l + 20)(l-2^/;) ^2^^2 



- [1 + 20 - (1 - 2?/;)cjV2 



d^ 



(1 



[1 + 20 
2?/;)wr2 



(1 



2^/>)a;2r2] 

2 



[1 + 20 - (1 - 2?/;)w2^2- 



■dp 



+ (1 - 2'ip)2r^LUUJ^rtdrdt - (1 - 2ip)2r^u^rtdrdp . 



(21) 



Unfortunately, this metric exhibits direct time dependence in some coefficients. The 
metric ()2ip is only reasonable in a strictly local sense. To preserve stationarity we 
approximate (|2ip by allowing only for small time intervals or equivalently for small 
angles of rotation. Then the differentially rotated PN metric (PT|) takes exactly the 
same form as only with u; = u{r) , and is thus also, to leading order, compatible 
with the Weyl gauge of the LP metric for the simplest PN case (j> = tp . 

Finally, let us note the result for the PN spacetime rotated fully differentially vi£@ 
Q = r2(r, z) which is a lengthy expression we write down in component notation: 

gtt = -(1 + 20) + (1 - 2i:)fl\^ , gtr = (1 - 2^)tr^nn^rt , fft^ = -(1 - 2i:yfl 
gt^ = {l-2^P)r^fin^,t, grr = (l - 2^j){l + r^n^.t^) , .g,^ = -(1 - 2^)r2rj,,i , 



g„ = (1 - 2tpyn,rn^,t' , 

g,, = il-2i')il + r^niy) 



g^^ = (1 - 2^y 



-{i-2^jyn,,t, 



(22) 



''Instead of lu we use the notation with Q to indicate that this can be understood, more generally 
than angular velocity, as just another allowed metric function. 
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Physically relevant however - again we need to meet the constraint of stationarity -, 
is (1211) taken as valid in a strictly local sense, that is 

ds2 = [-{l + 2(j>) + r^n^{l-22P)]dt^ + {1^2i^){dr^ + dz^+r^dip^-2r^nd<fidt). (23) 

After some algebra one sees that this is consistent with the Weyl form ([6]) and can 
be used as the starting point for a galaxy model in (Post-)Newtonian language. The 
analysis of the resulting dynamical equations via a 3+1 split is the final step and is 
going to be published elsewhere. 

3. DISCUSSION 

Motivated by the GR view on the galactic rotation curve problem, we have reviewed 
the outline of axisymmetric and stationary spacetimes as potential galaxy models. We 
have seen that a certain proposed model, the CT model, does not employ the most 
general axisymmetric and stationary spacetime form, the LP form, but belongs to a 
subclass thereof, the Weyl gauge. It is open in how far this restriction be connected 
with the difficulties of the CT model. 

If there really exist extra terms in a certain general relativistic approach then even- 
tually we should be able to pin down the differences by comparison to well-known New- 
tonian physics. As a first step we approached the CT metric in a (Post-)Newtonian 
language. We have shown explicitly that the classical Newton metric after (differen- 
tial) rotation cannot be brought to Weyl gauge, but is of a more general type. On 
the other hand, the result of (differential) rotation of the Post-Newtonian spacetime 
is consistent with the Weyl form in the case of (j) — tp . Thus the GR galaxy model 
considered requires a Post-Newtonian approximation. In a forthcoming publication 
we are investigating the consequences of these PN terms on galactic dynamics. 

This work was supported by the DFG under grants GRK 881 and GRK 1147. 
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